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Abstract. We consider the energy averaged two-point correlator of spectral 
determinants and calculate contributions beyond the diagonal approximation using 
semiclassical methods. Evaluating the contributions originating from pseudo-orbit 
correlations in the same way as in [S. Heusler el al. 2007 Phys. Rev. Lett. 98, 
044103] we find a discrepancy between the semiclassical and the random matrix theory 
result. A complementary analysis based on a field-theoretical approach shows that the 
additional terms occurring in semiclassics are cancelled in field theory by so-called 
curvature effects. We give the semiclassical interpretation of the curvature effects in 
terms of contributions from multiple transversals of periodic orbits around shorter 
periodic orbits and discuss the consistency of our results with previous approaches. 



1. Introduction 

It is currently accepted that progress on the long way to rigorous grounds of the 
Bohigas-Giannoni-Schmit conjecture pQ, stating the emergence of universal fluctuations 
in the quantum spectra of systems with classical chaotic dynamics, depends on the 
understanding of correlations among the actions of classical periodic orbits. This idea 
was pioneered in Berry's seminal paper [2], showing how to incorporate the most basic 
kind of action correlations into the semiclassical calculation of the spectral form factor. 
Berry's calculation, referred to as diagonal approximation, considers the contribution 
from orbit pairs with exactly the same action, i.e. orbits paired with themselves or 
their time-reversed counterpart in the case of time reversal symmetry. This calculation 
reproduces the leading order in a short-time expansion of the form factor obtained from 
random matrix theory (RMT) (or the leading order in an expansion in the inverse energy 
difference of the RMT spectral correlator). 

In seeking classical correlations responsible for the higher order terms Sieber and 
Richter succeeded [3] , invoking a subtle shift of the point of view: instead of looking for 
pairs of classical orbits with correlated actions, one constructs them for a given value of 
the effective Planck's constant. The mechanism underlying the formation of correlated 
orbit pairs is the exponentially close approach of classical orbits and their time reversed 
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partners everywhere except in the vicinity of a self crossing, as shown in Fig. [TJ This 
approach, giving the next to leading order contribution to the RMT form factor in the 
orthogonal case, was then generalized, leading to the semiclassical calculation of the 
RMT form factor for times smaller than the Heisenberg time T# [I]. Based on the 
analogy between semiclassics and field theory the form factor for times larger than T# 
could be obtained semiclassically in [5], although the underlying correlations remained 
unclear. 




Figure 1. Sketch of the orbit correlation mechanism described in [3]. 



The approaches dealing with times smaller than the Heisenberg time started from 
the density of states d(E), defined in terms of the spectrum {E n } as 

oo 

d(E) = Y,HE-E n ), (1) 

n=l 

which is approximated semiclassically by the Gutzwiller trace formula [6], 

d G (E) = d(E) + \^y"T^ im D p (E)e iSp{E)/h . (2) 

p 

It contains a smooth part, the mean level density d(E), and an oscillatory part expressed 
as a sum over classical periodic orbits p possessing the classical actions S p , the primitive 
periods T^ rim and the stability amplitudes D P (E); for their exact form see [6]. The 
Gutzwiller trace formula incorporates however the drawback that it is divergent, since 
it contains an infinite sum over an exponentially proliferating number of periodic orbits. 

This problem of the Gutzwiller trace formula renewed the interest in the 
semiclassical theory of spectral determinants. The quantum spectral determinant is 
defined as 

oo 

A(E) = Y[Sl(E,E n )(E-E n ), (3) 

n=l 

where Q(E, E n ) is a regularization factor without real zeros. The density of states, 
considered before, is proportional to the imaginary part of the logarithmic derivative of 
this determinant. The crudest semiclassical approximation for the spectral determinant 
is its representation as an infinite product over the classical orbits p occurring in the 
Gutzwiller trace formula. Equivalently, A(E) can be written as an infinite sum over 
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composite orbits, so called pseudo-orbits A with overall actions and stabilities given, 
respectively, by the sum and the product of the actions and stabilities of all possible 
subsets of periodic orbits with tia elements [7], 

A G (E) = B(E)e~ l ^ E) Y[exp{-D p (E)e iS *W /h } 

p 

= B(E)e- lnR{E) ^(-l) nA F A (E)e lSA{E)/h . (4) 

A 

Here, N(E) is the mean number of states up to energy E and B(E) a regularization 
factor. This form is divergent, as it is only based on the Gutzwiller trace formula. 
However, by directly imposing that the spectral determinant is real for real energies, a 
semiclassical theory of the spectral determinant can be constructed. Such an approach 
was developed by Keating and Berry [8] and obtained by Georgeot and Prange using 
Fredholm theory [9]. Thereby an improved (resurgent) semiclassical approximation to 
the quantum spectral determinant is obtained. It is given by 

A* (E) = 2UA% /2 (E) , (5) 

where A g h j 2 (-^0 * s obtained from A G (E) by including only pseudo-orbits with periods 
smaller than Th/2. Taking the real part assures that A Rl (E) is real for real energies. 
Moreover, truncating the sum over pseudo-orbits makes A Rl (E) convergent. In this 
paper we will also consider 

A Ra (E) = 29fJA G (E) , (6) 

imposing only the spectral determinant to be semiclassically real. 

Starting from the semiclassical approximation A G (E), Ref. [5] considered a 
generating function containing in the numerator and in the denominator two spectral 
determinants at four different energies. By taking derivatives with respect to two 
energies, identifying afterwards the two energies in the numerator with the two energies 
in the denominator in two different ways and adding the results of these identifications, 
they obtain the full RMT form factor. The idea then pursued in Ref. [10] was to use the 
representation (jSJ) for the spectral determinant instead of A G (E) and to consider only 
one identification of the energies. Due to complications from the finite cutoff of the sum 
defining A Rl (E), the authors of Ref. [10] used thereby A R2 (E). 

In this paper we consider the question, how to obtain semiclassically the RMT 
results for the energy averaged two-point correlator of spectral determinants, 

{ A(E + e/M)A(E- e/ 2l )}E 
' JM-oo <A(£) 2 ) E ^ ' 

starting from the semiclassical representation A R2 (E) and using the knowledge about 
action correlations between orbits and pseudo-orbits. The index (3 in Eq. ([7]) refers to 
the universality class, i.e. (3 = 1,2 for the orthogonal and unitary case, respectively. 
This problem was first addressed by Kettemann et al. [11] in the context of quantum 
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maps using diagonal approximation. The quantity Cp = i{e) was then reconsidered in 
diagonal approximation in Ref . [10] for continuous flows applying the methods of [5] . In 
both cases the calculations show perfect agreement with the universal RMT results for 
the unitary case 



Here the proportionality sign indicates that we focus on the universal part of Cp{e) and 
do not consider the semiclassical calculation of the prefactor. 

In the case of time reversal symmetry, however, Ref. [11] could not obtain 
semiclassically the RMT result 



In this paper we evaluate Cp(e) beyond the diagonal approximation. Using 
semiclassical techniques to account for action correlations and its natural extension 
to pseudo-orbit correlations we first demonstrate in Sec. 2 that the pseudo-orbit 
correlations, shown to cancel in the calculation of the spectral density correlator 
[5], do not vanish when one considers the simplest non-rational function of spectral 
determinants, i.e. Eq. ([7j). We show furthermore by explicit implementation of A Rl (E) 
that the problem, i.e. the apparent disagreement with RMT, is not solved with the 
correct truncation of the sum over pseudo-orbits. The third section is devoted to a 
corresponding field theoretical analysis of Cp(e), showing that the additional terms we 
find in our semiclassical analysis are cancelled in field theory by so-called curvature 
contributions. In the fourth section we then provide a new type of correlations between 
pseudo-orbits that allows us to semiclassically recover the RMT result in the unitary 
and, for the first time, also in the orthogonal case. After a discussion of the consistency 
of our results with previous results in the fifth section we finally conclude. 

2. Semiclassical calculation 

We semiclassically evaluate the correlator Cp(e) by means of A R2 (E). Upon substituting 
Eqs. (J3J) and (jSJ) into the numerator of Eq. (J7|), this leads to 



where ta,b = Ta,b/Th are the periods of the pseudo-orbits A,B in units of the 
Heisenberg time Tjj = 27fhd, and where we neglected highly oscillatory terms in E. 
We also expanded the actions of the pseudo-orbits around E. In the present context 
the application of the diagonal approximation is equivalent to taking the actions as 




(9) 




(10) 
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uncorrelated random variables, and therefore selecting A = B. We then can write the 
sum over pseudo-orbits as an exponentiated sum over orbits a. The Hannay and Ozorio 
de Almeida sum rule [12] is afterwards used to replace the sum over orbits, Wo? (■)> 
by where the minimal period T will produce an e-independent prefactor that 

cancels due to the denominator in Eq. (JTj). We finally obtain in diagonal approximation 
the first term in the squared brackets in 



[lire 



1 2/P 



l + C B (e) 



The function Cp(e) contains all universal effects beyond the diagonal approximation, 
consisting of both orbit and pseudo-orbit correlations. In the following we compute 
Cp(e) to leading order in 1/e. 

In the unitary case, the leading order in 1/e contribution to Cp = i{e) comes from 
the first pseudo-orbit correlations, namely, the "eight-shaped" diagrams, see Fig. [2] [5 J . 
To be more precise, we include in A a pseudo-orbit of order n and in B a pseudo-orbit 




Figure 2. Sketch of an "eight" diagram, an example of correlated pseudo-orbits [5]. 



of order n + 1 and vice versa. All orbits contained in A and B are assumed to be equal 
except for the ones forming an "eight" -orbit. This orbit consists of two components: one 
long orbit encountering itself at one point and two independent orbits encountering each 
other near the encounter of the first component. Both components are exponentially 
close to each other up to deviations in the encounter region. For the evaluation of this 
contribution we follow Ref. [5]: AS = su is the action difference of the two components 
with s, u being the difference of the stable and unstable coordinate in a Poincare surface 
of section placed in the encounter-region of duration t enc — •jj In J with a classical 
constant c. Furthermore the weight function counting how often the long orbit of length 
T has an encounter, is given by wt (s, u) = ^^wfy > where Vt(E) is the volume of the 
energy shell at energy E. Employing these expressions we obtain 

Cjgf(e)= -2 f C dsdu ^ ^ WT ( s ,u)e^ h > u e 2meT ^ 

J c J 2£ cnc 
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nie 

where the factor 2 in the first line results from the possibility to include the two short 
orbits in the A- or S-sum in Eq. (HDD- Using CjSf (e) in Eq. flJT]) gives a result which 
is inconsistent with the RMT prediction (JHJ). 

In the orthogonal case we must consider also the contribution C'pL 1 {e) from the 
"oo-shaped" orbit pair sketched in Fig. [JJ, apart from the "eight-shaped" contribution. 
Furthermore, in the latter the possible time-reversal operation leaving invariant each of 
the orbits must be taken into account. The calculations follow the same lines as for the 
unitary case, and we obtain 

Cfi\\e) =-—, (13) 
H Z7re 

C^(e) = — • (14) 
H me 

Again this result contains extra terms resulting from the "eight-shaped" diagram when 
compared with the corresponding RMT result (J9l). We note that the contribution 
Crtli(e) alone yields a result consistent with the RMT prediction. 

We can conclude that, although the analysis based on the semiclassical 
representation A R2 (E) of the spectral determinant produces the same structure as 
the RMT results (both are naturally ordered in powers in 1/e), the terms (fT2l [TBI) 
following from semiclassical loop contributions are not consistent with RMT and spoil 
the agreement found between RMT and the results in pJJ H] based on the diagonal 
approximation. Moreover, using the same methods we also obtain non vanishing 
contributions of higher order in 1/e [TBI fT4]. 

A natural candidate to solve the discrepancy is to use A Rl (E), Eq. (jSJ), in the 
semiclassical calculation instead of A R2 (E). The required truncation of the sum over 
pseudo-orbits can be exactly incorporated by means of the representation 

i r T H/2 poo 

(Tjj/2 — Ta) = — / dx dke ik{x - TA) (15) 
2tt Jo > /_ 00 

for the step function |15j . Including Eq. (fT5l) in the semiclassical calculation presented 
above, we end up with an integral which must be solved numerically [13]. The result 
clearly shows that the extra contributions coming from the "eight- shaped" diagram will 
not be corrected by using the fully resurgent expression (jSJ) for the semiclassical spectral 
determinant. Furthermore we rule out that the additional terms result from replacing 
in A Kl (E) the smooth step function 

1 f°° 

Q a (Ttf/2 -T A ) = -= dxe~ x2 , (16) 

V 71 " J(T A ~T H /2)/a 

suggested in [8], by a sharp cut off. More explicitly, using a Sommerfeld expansion 
we showed that the smooth cut off only adds to the result with sharp cut off terms 
proportional to o n with n e IN. These additional terms vanish when a is sent to zero 
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at the end of the calculation, meaning that the two limits cr — > and the semiclassical 
limit H — > commute. 

Obviously, the problem requires a better insight into the meaning of the pseudo- 
orbit correlations. Such an insight is gained through the comparison with a 
corresponding field-theoretical approach. 



3. Field-Theoretical calculation: the curvature contribution 



In order to better understand the difference between the semiclassical and the RMT 
results, we calculate Cp(e) by field-theoretical methods [16]. Our approach consists 
of transforming the non-perturbative field-theoretical expression for Cp{e) into a 
perturbative expansion. As we will see, the difference between the semiclassical and 
the RMT results has a direct correspondence in the field theory. 

We start with the non-perturbative evaluation of the integral following from field- 
theoretical considerations, which is given by [1] 



Ctfejoc-l f d^S)e~^ s , (17) 

with S G Sp{4)/ (Sp(2) * Sp(2)) in the orthogonal and S G U(2)/ (17(1) * 17(1)) in the 
unitary case. Parameterizing the latter space in terms of angles 9 and 0, we obtain in 
the unitary case 

C^ =2 (e)oc— I dcj) I dcos9e- t7Tecos9 



An _ 
sin 7re 



7re 

In a similar way, we obtain in the orthogonal case 

cosTre sinvre 

C/3=i(e) oc - —g-, (19) 

[Tie) {ire) 

in agreement with C^ MT (e). A perturbative evaluation of the field-theoretical 
expressions, in the form of a power-series in 1/e that is comparable to semiclassics, is 
obtained by replacing the integration variables on the sphere by stereographic projection 
variables in the complex plane. We therefore rewrite S as 

s I it l ] 7 1 (20) 




-I2//3 



with 




12//3 

and 7? G C in the unitary and 

B = f 1 1 ~~ 2 I (22) 
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with Bi, B 2 G C in the orthogonal case. We used here the dagger to indicate the 
Hermitian conjugate of a matrix and the power —1 for its inverse. 

The effect of this parameterization is twofold: On the one hand this variable 
transformation has a non-trivial Jacobian corresponding to a curved measure, 

4 [dB] 



d»(S) = - (23) 



with B = \B\, [dB] = dBdB* in the unitary and B = y|5i| 2 + \B 2 \ 2 , [dB] = 
dB\dB 2 dB\dB 2 in the orthogonal case. On the other hand it leads to a particular 
form of the phase in the exponential 

TrS= |^l + 2^(-l) fc £ 2fc j . (24) 

In this way and for this particular parameterization of the supermanifold where 
the field-theoretical calculation is defined, we can unambiguously identify two kinds 
of terms in the 1/e expansion of the correlator obtained by using in the unitary case 
J °° dB \B\ 2k Q- 2me \ B \ oc (l/e) fe+1 and the corresponding relation for B\ and B 2 in the 
orthogonal case. We will call these terms simply phase contributions if the factor in 
front of the exponential in the last integral results from expanding the right hand side 
of Eq. ff24l) and curvature contributions if it results from expanding the right hand side 

of Eq. (BHD i n P owers °f 

The key point is that we find agreement between the field-theoretical and the 
semiclassical results (112W14p . by evaluating Eq. ( fT7l) and ignoring the terms originating 
from the curvature of the measure, i.e. replacing (1 + B 2 )~ 4 ^ in Eq. ( 1231) by a constant 
[T3~t [14] . Clearly, what we are missing in the semiclassical calculation are the analogs 
of the terms that can be derived from the curvature. This means that a new kind of 
pseudo-orbit correlation is lacking which corresponds to the curvature effects. Moreover, 
one has to show that this new correlation is consistent with existing semiclassical results 
showing agreement with RMT. 

It is important to notice that curvature effects did not appear in the case of the 
spectral generating function in Ref. [5] because in the corresponding field-theoretical 
analysis the curvature factor had to be taken there to the r-th power with the Replica- 
index r and was thus vanishing in the Replica-limit r — > 0. In practice, this observation 
can be extended to any measure of the fluctuations which is expressed as an homogeneous 
rational function of spectral determinants. 



4. Semiclassical interpretation of curvature effects 

In this section we present the missing pseudo-orbit correlations which provide the 
curvature effects in the integration over the curved manifold in the field-theoretical 
approach and thereby consistency with RMT. The essential step is to lift the constraint 
that one of the short periodic orbits does not surround the other one more than once. 



The semiclassical origin of curvature effects in universal spectral statistics 



9 



This assumption is implicit in the calculation of the "eight-shaped" diagrams in Ref. [5] 
when one assumes that both short orbits have periods larger than the encounter time. 

To be precise, we consider as components of a pseudo-orbit two periodic orbits 
(dashed and full line in Fig. [3]) which include several transversals around a common 
short periodic orbit and differ in the number of repetitions by one. We then pair the 
short periodic orbit and the orbit surrounding it (k — 1) times with the orbit surrounding 
it k times. This situation occurs in the case, when the encounter time t enc is longer than 
the shortest of the periodic orbits, i.e. when the encounter overlaps. 




Figure 3. Sketch of a pair of pseudo-orbits A, B correlated in a way that one of 
the orbits (dashed-dotted) has a period smaller than the encounter time. The two 
long orbits (dashed and full line) differ in the number of transversals around the short 
periodic orbit by one. 



For the evaluation of this contribution, we again need the action difference between 
the two orbits, obtained as AS = su [17J up to a negligible correction and the weight 
function as a function of s, u that is given by wtt„ (s, u) = t^te& — expressed through 
the periods of the two short periodic orbits, T and T p . Using these results and applying 
the sum rule with respect to T and T p , we find for Cp=2(e) the additional "curvature" 
contribution 

/C POO irp pt cnc irp 

dsdu ^ ^ e m)su e 2 m e { T + T p )/T HwTTpM 

c Jt cnc * JO 

= — , (25) 
me 

thus precisely canceling the undesired term (I12p coming from the "eight-shaped" 
diagram, and allowing for the full explanation of the RMT result (|HJ) in purely 
semiclassical terms. 

In the orthogonal case, one obtains in a similar way 



Cfll{e) = — , 
tire 



(26) 
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canceling again the extra term (TT5j) and yielding the universal RMT result (151) . 

We could thus identify the semiclassical analog of the field theoretical curvature 
contributions. As in the field-theoretical case the splitting between the two diagrams 
analyzed in the semiclassical calculation in Sections 2 and 4 is completely arbitrary; 
the final result, the sum of the two is however always the same. This can be taken 
into account by considering as reference orbits in the calculation of the weight function 
always the two short orbits as we did above Eq. (125]) and by allowing for an arbitrary 
length of one of the two periodic orbits as long as the other one is longer than the 
duration of the encounter. This yields 

Cjg'(e) + C—{e) = - 2 f dsdu jf ^ £ ^m™^*^^ (S) u) 

dsdu J ±LJ ^e^e™< T+T *y T »w T , Tp (s,u) 

= 0, (27) 

i.e. the sum of the two contributions calculated in the sections 2 and 4. The splitting 
of the two contributions in our semiclassical calculation showed however the one- 
to-one correspondence to the corresponding field theoretical contributions after the 
stereographic projection. 

We want to emphasize here the different behavior of an "eight" orbit and the "oo" 
orbit pair in Fig. [TJ for an "eight" orbit with overlapping encounters, one of the short 
orbits winds around the other and thus yields another nontrivial contribution, that is not 
yet contained in the diagonal terms. In the case of an "oo" orbit pair with overlapping 
encounters, however, the orbit and its partner transverse the whole orbit in the same 
direction, which is already contained in the diagonal contribution. 

By this calculation we demonstrated that overlapping encounter regions, never 
considered before in evaluations of energy-averaged quantities neglecting the Ehrenfest- 
time dependence |18j . can lead to non- vanishing contributions even for vanishing 
Ehrenfest-time. One question, however, has to be addressed at this point: When are 
effects of encounter overlap important and when can they be neglected? This point will 
be answered in the following section. 

5. Consistency with former results 

Also in former semiclassical calculations to reproduce RMT-results diagrams involving 
orbits surrounding a short periodic orbit occurred. Here we discuss the consistency of 
our results with those approaches. 

5.1. Approaches involving orbits (but no pseudo-orbits) 

In this case an overlap of encounters can only happen for different encounters, that 
means, it is not possible that the ends of one encounter come close to each other, 
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Figure 4. Encounter overlap for different encounters 

because otherwise there exist only disconnected partner orbits. One simple example for 
a possible diagram, where a longer orbit surrounds a short periodic orbit is shown in 
the Fig. H 

In this case, however, a calculation similar to that for the spectral form factor for 
finite Ehrenfest-time [H] shows that the resulting contribution contains, for an arbitrary 
time of the surrounded periodic orbit T p , just the configurations of two independent 
encounters of two orbital parts and one encounter of three parts. Such configurations 
were already taken into account in [1]. In particular, this means that we do not find 
a different result by also taking into account a configuration where a periodic orbit is 
surrounded by a longer one many times. 

5.2. Approaches involving pseudo-orbits 

In Ref. [5], the contribution from "eight-shaped" orbits vanishes also, but for a different 
reason: the contributions arising from different possibilities to include the "eight- 
shaped" orbit in the pseudo-orbit sums cancel (giving agreement with RMT). However, 
this method is not applicable in general, as was shown in this paper. Still this will not 
pose a problem for measures of the spectral fluctuations based on the density of states, 
since they involve homogeneous rational functions of spectral determinants. 

6. Conclusions 

We have shown that the semiclassical approach from Ref. [5] applied to the two-point 
correlator of spectral determinants (and we believe this also holds for any measure 
involving only products of spectral determinants) produces spurious terms coming 
from pseudo-orbit correlations, spoiling agreement with the universal RMT results. 
We argued that the appearance of extra terms in the correlator is not solved by the 
strict use of the resurgent spectral determinant, and that a new additional mechanism 
is required. This contribution consists of pseudo-orbit correlations represented by 
encounter diagrams with overlapping encounters. We have shown that these effects 
are the semiclassical analogue of the curvature contributions in the field theoretical 
approach, thus providing a semiclassical interpretation of the latter. 
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Our results can be extended in various directions: First, concerning the study of 
action correlations between orbits it would be interesting to further analyze the new 
diagram type with an overlapping encounter identified in this paper: Are there other 
quantities apart from the product of two spectral determinants, where this correlation 
is essential to reproduce the RMT-results or do similar correlations lead to further 
undiscovered contributions? Second, as the correlator of two spectral determinants can 
be used to describe certain features of localization [19], it would be interesting to use 
the field theoretical calculation presented in Ref. [19] as a guideline for a semiclassical 
approach to localization. 
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